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I. INTRODUCTION 

For homogeneous triaxial ellipsoid with semiaxes A, B, C and density p the next relations are valid (Landau and 
Lifshits 1975).The gravitational potential at the inner point {X, Y, Z ), -A < X < A, -B < Y < B, -C < Z < C, 
is 

UMX,Y,Z) =npGABC fl _ ^ _ ^\As. 



Qs = ^{A^ + s){B^ + s){C^ + s). 



^> ' The potential energy of the homogeneous triaxial ellipsoid is: 



is 



Here Meu = A/3tt p AB C is the ellipsoid's mass and G stands for Newtonian constant of gravitation. Note that 
, gravitational energy of self-gravitating body is of negative sign but we loosely write all Ws with positive sign. In 
general case A ^ B ^ C the integrals in (^ and (||) are expressed only in terms of the incomplete elliptic integrals 
^ ^ and this precludes any detailed analysis. However, if we consider only potential at the center of ellipsoid, Ueiiifi, 0, 0), 
I > then we get the remarkable relation: 

o , 

+5 ■ WUM = -, (3) 

J> ■ valid for any values of semi-axes. We shortly refer to this relation (|^) as WUAI-iatio. 

' Recently Seidov and Skvirsky (2000a) presented the gravitational potential and the potential energy for the homo- 
. geneous rectangular parallelepiped (hereafter RP) which allows to analyse this WUM-ratio for the new class of the 
homogeneous self-gravitating bodies. 

In this paper, we show in sections II and HI that the value of WUM for RPs has minimal value .395437 (see Eq. 
(p^)) for the cube (all three dimensions of RP equal to each other), tends to 1/2 for one dimension of RP far larger 
larger than two others (long thin "stick" with square cross-section), and tends to ^ — g ~ -421673 as one 

dimension is far less than two others (thin square "plate"), see Fig. 1. 



Also, in section IV we discuss the WUM-idXio for homogeneous gravitating bodies studied recently by Kondrat'ev 
and Antonov (1993). We show that values of VFt/M-ratio for the homogeneous summetrical lenses are in the interval 
from 128/105 TT, for infinitesimally thin symmetrical lens, to 17/20, for two homogeneous equal spheres just touching 
each other. 



In section VI we analyse WU M-idXio for spherical polytropic stars, and show that VF[/M-ratio varies from 2/5 to 



3/32 TT for polytropic index n varying from to 5. 



In the section |VH| we discuss the interesting class of two-phase spheres and show that unlike the polytropes, in this 
case the WU M-ratio's interval is larger: it is possible to get very small values of WU M if the ratio of two densities 
(1 = Pil P\ is large enough and if the relative value of core's radius is not too small. 



At last, in the sections VIII and tX we consider another two simple classes of non-homogeneous bodies both allowing 
analytical treatment. 
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II. POTENTIAL AT THE CENTER OF RP 

Using results by Seidov and Skvirsky (2000a) we write down the gravitational potential at the center of the homo- 
geneous RP with density p and with edge lengths 2a, 26, 2c: 



C/flp(0,0,0)-4Gp[a61n^ + 6cln^ + crflng- 
arctan —, —b"^ arctan f-§ — arctan ^1. 

ad b a c a' 



(4) 



Here d = (a^ + 6^ + c^)^/^ is the main diagonal of RP. 
Three particular cases are of the larger interest: 
a) cube corresponding to case c — b — a, 



241ni±^-2. 
V2 



9.52017Gpa2; (5) 



h)long thin stick with square cross-section corresponding to case a >> b = c: 

Usuck (0, 0, 0) = G p fe2 (_8 In ^ + 12 _ 2 ^ + 4 In 2); (6) 

a 

c) thin square plate corresponding to the case a << b = c: 

Upiate{0,0,0) - Gp[16a6 ln(\/2 + 1) -2T:a% (7) 

III. POTENTIAL ENERGY OF RP 

According to Seidov and Skvirsky (2000a) the gravitational potential energy of the homogeneous rectangular par- 
allelepiped is equal to: 

WRP ^Gp^[f{a, b, c) + /(6, c, a) + /(c, a, b)]-J{a, b, c) ^ c^a^ + 040^ + c^a^ + csa^; 



C5 = f|; C4 = i(d-dl-rf3)-i|^ln(i^i^l^-i|^ln(i^5^^ 

C3 ^ -64^ arctan^; d = Va^ + 6^ + c2; dl = Va^ + b^; d3 = Va^ + c2; 

C2 = '^{dl-d) + ^{d'S-d) ~ 16 6c2 Inf^ - 16 6^0 In 



(8) 



A. Potential energy and WUM-ratio of cube 

From (^, taking c = 5 = a, we get the potential energy of homogeneous cube with edge length 2a: 

= 32GpV{ 2^^/^-i +f+ln[(V2 -1) (2- V3)]} = 30.117 Gp^a^. (9) 

From (^) and (|^) we get T^JJM-ratio for homogeneous cube: 

^^^^l^ + ^+ln[(V2-l)(2- V3)] , , 

WUM^ubc = 2 § 2 \11 11 ^ = .395437. (10) 

24 In ii# - 27r 

v2 
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B. Potential energy and WUM-ratio of thin long stick 



We take a >> b = c that corresponds to the case of the thm long stick with the square cross-section. Leading term 
in expansion of WRP (||) gives the potential energy of the thin long stick: 



tick 



|Gp^a6^1n^. 



From this and (m we get for stick: 



WUM = 



stick 



8 a b^UsUckiO, 0,0) 



(11) 



(12) 



C. Potential energy and WUM-ratio of thin square plate 



Taking one of RP's dimension infinitesimally small, a — + 0, we get, from Eq.(^, the potential energy of the thin 
rectangular plate. If we additionally take b = c, then we get the potential energy of the thin square plate (a << b): 



W^i = 64 



Gp^b^a^ ^ln(V2 + 1) - j = 47.5714 G p^&^a^. 



From (0) and (|T|) we have another limit for value of WUM: 



WS 



%/2- 1 



8a52C/p,(0,0,0) 2 6 1n(\/2 + l) 



.421673. 



(13) 



(14) 
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FIG. 1. WUM for RP with square cross section 2 6x26 and length 2a. Abscissas are values of h/a and ordinates are values 
of relation WUM — W/U{0) M that is ratio of gravitational potential energy of RP to product of gravitational potential at 
the center of RP by mass of homogeneous RP. WUM has minimum for cube {b/a = 1), tends to 1/2 at b/a — > (thin long 
stick) and to i — g 1^(^2+1) ~ -421673 (dash line ) at b/a 00 (thin square plate). For homogeneous ellipsoid, WUM = 2/5, 
solid line. 



D. Relation between potential energy, gravitational potential and mass of RP 

General behavior of relation between the potential energy, the gravitational potential at the center, and mass of 
the homogeneous RP with two equal edge- lengths is shown in the Fig. 1, which is a result of numerical calculation 
by the formulas (^) and (^. For homogeneous ellipsoid WUM = 2/5, see (^, solid line in Fig. 1. 



IV. HOMOGENEOUS SYMMETRIC LENSES 



Recently Kondrat'ev and Antonov (1993) (hereafter KA) have obtained the analytical formulas for the gravitational 
potential and the gravitational energy of some axial-symmetric figures, namely homogeneous lenses with spherical 
surfaces of different radii. In forthcoming paper (Seidov and Skvirsky 2000b) we present some new solutions for 
homogeneous bodies of revolution. 

Here we present the review of WUM for most suitable kind of those bodies, discussed by KA, namely the homogeneous 
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symmetrical lenses. A segment of sphere, or a planoconvex lens, is obtained by cutting a sphere with a plane. If R is 
a radius of a sphere, h is height of segment, and 2 a is a radius of segment's base, then a = h R — h? . A symmetric 
homogeneous lens (SL) is obtained by placing together the bases of two identical segments of sphere. 
According to KA we have for the gravitational potential at the center of such SL: 

UsLiO)^^^{hR + R^-- + -^—j^]. (15) 



The gravitational energy of SL is: 



WsL 



~ 9 



lOR^a - fR^a-^ - ^a^ + nh^{fh~2R- ^) + 2 i?M - 6 (i? - /i) arctan 



And the total mass of the homogeneous SL is: 



(16) 



MsL^^^{iR~h). (17) 
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FIG. 2. WUM for SL. Abscissas are values of h/R, ratio of height of half-lens to radius of sphere and ordinates are values 
of relation WUM = W/U{0)M. At h/R (thin spherical lense), WUM ^ = .388035. At h/R 1, (full sphere) 
WUM ^2/5. 



We may consider WUM-ratio for such bodies as Wsl/Usl Msl- Defined so, WUM for the homogeneous symmetric 
lenses has general dependence on parameter h/R, according the formulas (|l^), (|l6|) and (p7|), as shown in the Fig. 
2. Note that in the limit h/R ^ we have infinitesimally thin symmetrical lens with radius of curvature tending to 
infinity, and we have: 

1 28 

Limit[WUMsL]h/R ->o = TTT^ = .388035. (18) 

' 1U5 TT 

Interestingly, this infinitesimally thin round lens docs not coincide at all with the case of the infinitesimally thin 
quadratic plane which has the much more larger value of WUM — .421673 (see Eq.(|4|)). In another limit h/R 1 
we have a full homogeneous sphere, and evidently: 

Ltmtt[WUMsL]h/R^>i = 2/5. (19) 

There is still another analytical case at h 2 R when we have two homogeneous spheres touching each other so 
that distance between centers of spheres is d = 2 i?; then "potential at the center of SL" is C/sl(0) = 2 potential 
energy is sum of two terms: proper potential energy of each spheres 2 | , and interaction energy, ^^^^ or 
we have: 

WUMsL\h=2R^{l + l)/2^^. 



V. HETEROGENEOUS SPHERICAL BODIES 



Now we consider the problem of WU M-ratio from another point of view. Abovementioned homogeneous gravitating 
bodies (ellipsoids, right parallelepipeds, and symmetrical spherical lenses) differ from each other only by their forms 
and so M^J7 M-ratio may be referred to as form-factor. As it is evident, M^t/ M-ratio should be also function of density 
distribution over the body. If we confine ourselves by spherically-symmetric distribution of density p(r), then we have 
general expressions for the gravitational potential at radius r: 
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Gm{r) 



J r 



central potential: 



potential energy: 



and mass: 



U{r) = —+ A-kG p{r)rdr; (20) 



C/(0) = / AiTGp{r)rdr; (21) 
Jo 



1 

W=- U{r) dm{ 
2 Jo 

An p{r)r'^ dr; M= / 4 7rp( 
Jo 



(22) 



dm = An p{r)r^ dr; M= / Anp{r)r^dr. (23) 

Jo 

Prom these expressions we write down WUM for spherical-summetrical heterogeneous bodies as: 



W 1 U ^ 1 



M 



WUM = —— = U=— U{m)dm. (24) 
Z7(0)M 2C/(0) M Jq ^ ' ^ ' 

As a result, T/KC/Af-ratio of spherically-symmetric bodies is reduced to the ratio of the mean value of nionotonic 
function, U{m), to its particular value, U{0), (ratio being additionally divided by 2). The boundary values of this 
ratio can be found pure mathematically for any given class of functions p (r) . We will not deal with this abstract 
(though interesting) problem; instead, in the next sections we consider two cases of more or less realistic bodies, 
namely polytropes and two-phase spheres. 

VI. POLYTROPES 

One case of heterogeneous bodies which apparently should be considered first is the case of classical polytropic 
stars. Using formulas from Chandrasekar's (1957) classical text we have (some of these formulas are valid not only 
for polytropic stars but we do not stop on these details) the next relations. 

a)central potential is expressed via other parameters of star as follows (Chan 100/85= Chandrasekhar (1957), p. 100, 
Eq. (85)): 

, . , .Pc GM 
[/p(0) = (n + l)^ + --. (25) 

Pc -K 

Here n is the polytropic index, Pc and pc are central values of pressure and density, M and R are the mass and radius 
of the star. 

The next formulas include the parameters of the Lane-Emden function (LEF) at the first zero point: 



Central pressure is (Chan 99/80,81): 



= ei'<0; ;xi = -^2^i'. (26) 



P - 1 (27) 

^~ 47r(n + l)(^i')2 ■ ^ ' 



Central density pc is related with the mean density p of the star as follows (Chan 78/99): 

1 ^ _ 

Combining these formulas we get the final expression for the central potential of polytropic star: 



(28) 



5 



Up{0) = 1 



1 \ GM 



1 oi 



R 



(29) 



h)potential energy 

We have for polyropic star the famous formula (Chan 101/90): 



c) WUM-ratio: 



_3_GM2 
R ' 



(30) 



WUM„ 



5-71 1 + ^ 
Ml 



(31) 




12 3 4 5 

FIG. 3. WUM for polytropes. Abscissas are values of polytropic index n and ordinates are values of the WU M-ra.tio. 

A. WUM-ratio for polytropes 

Using values of parameters of polytropic stars given in (Chan, Table 4, p. 96) we calculated the values of WU M- 
ratio for values of polytropic index n from to 4.9, see Fig. 3. However the limit n 5 and in general region of values 
of n close to 5 should be considered separately. First we note that atn^5, ^i-->oo and there is indeterminacy 
of ^ kind in formula for WUMp Eq. (|l]). To solve this problem we use results of Seidov and Kuzakhmedov (1979) 
who obtained, in particular, the dependence of i^i for n close to 5: 



(0 < (5-n)/5 « 1), a 



32^3 1 



TT 5 — n 



(32) 



Also it was shown by Seidov and Kuzakhmedov (1979) that for values of n close to 5, the parameter /ii is an increasing 
function of n: 



Ml 



\/3(l 



1 1 



12 5~n 



(33) 



that means that fii has minimum at n about 4.82. This may or may not lead to minimum of WUMp as function of 
n. At point n = 5 using LEF of index 5: 



n = 5, e{0 = {l + ^e)-'/', 



(34) 



we get 



n = 5, WUMp ^ —n^ .294524. 



(35) 
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Additionally wc have two other analytical results for WUM of polytropes : at n = 0, WUM = 2/5 and at n = 1. 
WUM = 3/8. We recalculated the parameters of polytropes in the interval of n from 4 to 5 and found no minimum 
for WUM, see Fig. 3. However we note that our boundary values differ from ones calculated by Jabbar (1993) in the 
sense that ours arc less than his. As one example, at n = 4.7 Jabbar gives = 54.810686 as zero of LEF, while our 
calculations using Mathematica's command NDSolve give 0(54.810686) = —4.5910" and our value of is in interval 
between ^ = 54.8098 (at this point 9 = 4.98 * 10-8) and ^ = 54.8099 (at this point 6* < 0). In general our values of 
are less than Jabbar's. 



VII. TWO-PHASE SPHERE 



If pi and P2 = qpi are densities in envelope and core of sphere, and R and r = xR are total radius of sphere and 

radius of core, then we have next relations: 
a) gravitational potential at the center: 

U2ph{0)=2 7rGpiR^[l + {q-l)x% (36) 

h)potential energy: 

W2,H = ^^^Pi' R'[l + l{q- 1) x' + {q- l){q - ^) x']; (37) 

c) total mass: 

M2ph = ^PiR''[l + {q-l)x% (38) 

d) WUM - ratio: 

WUM2,, = r]}x^::,^^iif . (39) 




0.2 0.4 0.6 0.8 1 

FIG. 4. WUM for two-phase spheres. Abscissas are values 
WUM-ra.tio. Dash line shows the locus of minimuma of WUM. 



of relative radius of core, and ordinates are values of the 



In Fig. 4 we present dependence of WU M-vatio, for two-phase spheres with various values of q, as function of 
relative radius of core x = r/R. Curves from uppper one to lower one correspond to values of g = 3, 5, 10, and 
20, respectively. Both values of minimuma of WUM and their "places", corresponding values of x, are decreasing 
functions of q. For values of x corresponding to minimuma of WU M the next simple equation is valid: 

6q'^x^ + 5q{l-x)x^{l + 2x)-{l-xf{4: + 7x + 4:x'^) = 0. (40) 



or 



(l-x) (A/96/a; + 97 + 28a; + 4a;2-5-10V^ 
12 a; 

We used these formulas to calculate the dash line in Fig. 4. 



9= ^ ^2 -■ (41) 
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VIII. STEPENARS 



Here we briefly consider the case of simple spherically-symmetric density distribution law allowing analytical ex- 
pression for WUM. We take: 



p{r)^p,ll~r/R)\ 



(42) 



where pc is central density and z/ is a free parameter. By historical reasons we refer to the gravitating bodies with 
density distribution ( ^2|) as "stepenars", (Seidov, Kasuniov, and Guseinov 1971). 
We have: 



mass: 



mean-to- central density ratio: 



central potential: 



potential energy: 



WUM-ratio: 



M = 



8 TT PrR^ 



(l + z/)(2 + i/)(3 + z/)' 

P_ ^ 6 

Pc (l + J/)(2 + z/)(3 + z/)' 

AnGpcR^ 



U{0) 



W 



(l + z/)(2 + z/)' 



(l + z/)2(2 +iyf{3 + 2iy){5 + 2v)'' 
(3 + i/)(8 + 5i/) 



WUM^ 



4(3-t-2i^)(5 + 2!y)' 



(43) 



(44) 



(45) 



(46) 



(47) 




FIG. 5. WUM for stepenars. Abscissas are values of parameter v in density distribution law p(r) = pc (1 — r/K)", and 
ordinates are values of Vl^i7M-ratio. Broken line shows asimptotic value of WUM 5/16 at u <x. 

In Fig. 5 the dependence of WUM-Tatio for stepenars as function of parameter z/ is shown. In spite of large variation 
of matter concentration from p = const at z^ = 0, to p/pc ^ at z^ — > oo, WUM-ia.tio again as in polytrope's case 
lies in rather limited interval from 2/5 to 5/16. 



IX. ALPHARS 

Here we consider another example of "exotic" but simple density distribution, Seidov and Seidova (1971): 

p{r)^Pc{r/Rr, 0<a<2. 

We have: 

2-a 



WUMc, 



5-2a 



(48) 
(49) 
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X. DISCUSSION 



There is a rather classic problem of looking for general theorems of stellar structure, see e.g. chapter 3 in the classic 
text Chandrasekhar (1957). The problem considered in this paper may be also referred to as that dealing with general 
structure of celestial self-gravitating bodies. 

We start from interesting observation on one constant ratio, namely, (potential energy W)/((central potential U ) x 
(total mass M)), in homogeneous ellipsoids and then try to look for behavior of this ratio for another homogeneous 
bodies: rectangular parallelelepipeds and symmetrical lenses. We found that in both cases WU M-ratios are confined 
in rather narrow interval. Suprisingly, dependence of WUM for homogeneous rectangular parallelepipeds (RP) on edge 
lengths ratio is non-monotonic: it has minimal value .395437 for cube while any deviation from cube form to prolate 
RP (one dimension being smaller than two others) or elongated RP (one dimension being larger that two others) 
leads to the increase of value of WUM. In this respect the behavior of homogeneous rectangular parallelepipeds is 
quite unlike the behavior of homogeneous ellipsoids and there is still some mystery even to authors. 
As to the homogeneous symmetrical lenses (SL) by Kondrat'ev and Antonov (1993), here the dependence of WUM 
on parameters of SL is monotonic, however in this case there is also some suprise in the sense that in the limiting 
case of thin symmetrical spherical lens the H^[/M-ratio's value, (128/105 tt), differs radically from the case of the 
infinitesimally thin quadratic plate with WUM = 1/2. 

Then we look for the non-homogeneous however spherically symmetric bodies and found that for the polytropes with 
polytropic index n in the interval — 5, WUM-Ta.tio again lies in narrow interval from 2/5 to 3/32 tt. However for 
two-phase sphere with large ratio of densities q = P2/P1 it is possible to get very small values of WUM. The physical 
reason of it is that if we put in the center of any spherical symmetric star a (very) small but dense spherical body 
then central potential may be very large while total potential energy of star, being integral value, increases not so 
drastically. The effect of the strong variation of density in the center of star, the "first-order phase-transition", is 
known since pioneer works of W.H. Ramsey (1950). 

In last two paragraphs of paper we consider pure mathematical toy models in further attempts to understand the 
behavior of the WU M-ratio. We conclude this discussion with notice that central-to-surface potential ratio U (0)/C/(i?) 
(among other "global" characteristics of the celestial self-gravitating configurations) is also worth studying. For 



polytropes, U{0)/U{R) = 1 +^i//xi, see section VI 
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